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ABSTRACT 

From observations collected with the ESPaDOnS and NARVAL spectropolarime- 
ters, we report the detection of Zeeman signatures on the classical T Tauri star (cTTS) 
BP Tau. Circular polarisation signatures in photospheric lines and in narrow emission 
lines tracing magnetospheric accretion are monitored throughout most of the rotation 
cycle of BP Tau at two different epochs in 2006. We observe that rotational modula- 
tion dominates the temporal variations of both unpolarised and circularly polarised 
spectral proxies tracing the photosphere and the footpoints of accretion funnels. 

From the complete data sets at each epoch, we reconstruct the large-scale mag- 
netic topology and the location of accretion spots at the surface of BP Tau using 
tomographic imaging. We find that the field of BP Tau involves a 1.2 kG dipole and 
1.6 kG octupole, both slightly tilted with respect to the rotation axis. Accretion spots 
coincide with the two main magnetic poles at high latitudes and overlap with dark 
photospheric spots; they cover about 2% of the stellar surface. The strong mainly- 
axisymmetric poloidal field of BP Tau is very reminiscent of magnetic topologies of 
fully-convective dwarfs. It suggests that magnetic fields of fully-convective cTTSs such 
as BP Tau are likely not fossil remants, but rather result from vigorous dynamo action 
operating within the bulk of their convective zones. 

Preliminary modelling suggests that the magnetosphere of BP Tau extends to 
distances of at least 4 R+ to ensure that accretion spots are located at high latitudes, 
and is not blown open close to the surface by a putative stellar wind. It apparently 
succeeds in coupling to the accretion disc as far out as the corotation radius, and could 
possibly explain the slow rotation of BP Tau. 

Key words: stars: magnetic fields - stars: accretion - stars: formation - stars: 
rotation - stars: individual: BP Tau - techniques: spectropolarimetry 
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1 INTRODUCTION 

T Tauri stars (TTSs) are young low-mass stars that have 
emerged from their natal molecular cloud core. Among 
them, classical TTSs (cTTSs) are those still surrounded by 
accretion discs. CTTSs host strong magnetic fields thought 
to be responsible for disrupting the central regions of their 
accretion discs and for channelling the disc material to- 
wards the stellar surface along discrete accretion funnels. 
This process is expected to play a key-role in setting the 
angular momentum evolution of Sun-like protostars (e.g. , 
Konigll Il99ll; IShu et all Il994l ; ICameron fc Campbelll 1 19931 ; 



Romania etal. l l 2004h as-^ell as their internal sMuTe" 



One of the prototypical and most extensively stud- 
ied cTTS is BP Tau in the Taurus star formation region. 
BP Tau hosts multi kG magnetic fields on its surface, de- 
tecte d both in photospheric lin es (using Zeeman broadening, 
e.g., IJohns-Krull et al.lll999ah and in emission lines forming 
at th e base of accretion funne l s (using Zeeman polarisation , 
e.g., IJohns-Krull et al.lll999b1 ; IValenti fc Johns-Krullll2004l) . 
However, the large-scale topology of the magnetic field is 
not well known, preventing detailed studies of how the star 
magnetically couples to its accretion disc and how efficiently 
angul ar momentum is transferred from th e disc to the star 
(e.g., iJardine et a"i1l2006l ; lLong et al.ll2007l ). 

Large-scale topologies of stellar magnetic fields can be 
constrained using Zeeman polarisation signatures in line 
profiles and the information they convey about the field 
orientation. By measuring the polarisation spectrum of a 
star and by monitoring its modulation throughout the whole 
rotation cycle, the large-scale topology of t he parent sur- 
face magnetic field can be retreived (e.g., iDonatil l200ll ; 
iDonati et aI1l2006bh . In the particular case of cTTSs, this 
requires Zeeman circular polarisation signatures to be de- 
tected, not only in emission lines - tracing the tiny areas at 
the footpoints of accretion funnels - but also in photospheric 
lines - forming over most of the stellar surface. 

This technique was recently applied to V2129 Oph, a 
bright, mildly accreting cTTS in the Ophiucus stellar for- 
mation region |Donati et al.1120071 hereafter D07). The mag- 
netic topology of V2129 Oph is found to be significantly 
more complex tha n the conventional dipole u sed in most 
theoretical studies l|Romanova et al ]|2003l. 120041 ); despite be- 
ing rather weak (0.35 kG), the dipole field component is 
nonetheless apparently capable of coupling to the accretion 
disc up to distances of 7 stellar radii. In particular, this re- 
sult brings further support to the original idea that magnetic 
coupling between the star and its accretion disc is able to 
control the rotation rate and the angular momentum content 
of the protostar. In this paper, we apply the same technique 
to BP Tau. 

Fr om a detailed spectroscopic study, IJohns-Krull et all 
( 1999a) conclude that BP Tau has a photospheric tempera- 
ture TLff of 4055± 112 K, a logarithmic gravity of 3.67±0.50 
and a logarithmic metallicity (with respect to the Sun) of 
0.18±0.11. Assuming a logarit hmic luminosity (with r espect 
to the Sun) of -0.03 ± 0.10 jGullbring et al.lll998l ). they 
derive that BP Tau has a radius of i?* = 1.95 ± 0.26 Rq. 
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The rotation period of BP Tau, as derived from photometric 
variability, is reported to be 7.6 ±0.1 d (|Vrba et al.lll986l : 
IShevchenko et all Il99ll ); discrepant estimates, including a 
sudden change of the light-curve period (from 7.6 to 6.1 d, 
ISimon et ai"1 [l990) are also reported, suggesting that the 
photometric brightness of BP Tau may sporadically include 
a non-stellar contribution, e.g., from the inner edge of the ac- 
cretion disc. Given the rotation period and the line-of-sight 
projected eq uatorial rotation velocity v sin i (~ 10 kms - 
according to IJohns-Krull et all Il999a. in good agreement 
with our own estimate of 9 kms -1 , see below), we find that 
R+ sini ~ 1.4 Rq and thus that i ~ 45°; 

Fitting the evolutionary models of lSiess et al.l l|2000l ) to 
these parameters, we infer that BP Tau is a 0.70 ± 0.15 Mq 
star with an age of about 1.5 Myr. Models indicate that 
BP Tau is still fully convective, making it thus significantly 
different from V2129 Oph (which recently started to build up 
a radiative core, D07) and thus especially interesting for our 
studjQ. Infrared excess and emission lines are conspicuous 
in the spectrum of BP Tau, indicating the presence of an 
accretion disc. Optical veiling, weakening the strength of 
photospheric lines by a significant fraction (typically 10% to 
50% depending on epoch and wavelength) also suggests that 
accretion hot spots are present at the surface of the star. The 
mass accretion rate of BP Tau is M ~ 3 x 10 -8 M Q yr -1 
(e.g., iGullbring et"ai"1l 19981 ). 

After presenting and describing our spectropolarimetric 
observations (Sees. [2] & [3]), we apply tomographic imaging 
on our time-resolved data sets and derive the large-scale 
topology of BP Tau at 2 different epochs (Sec. [4]) that we 
use to provide a first model of the magnetosphere and accre- 
tion funnels (Sec.0. We finally summarise our results and 
discuss their implications for our understanding of stellar 
formation (Sec. [BJ. 



2 OBSERVATIONS 

2.1 Spectropolarimetric data 

Spectropolarimetric observations of BP Tau were collected 
in 2006 at 2 different epochs, using mostly ESPaDOnS 
on the 3.6 m Canada-France-Hawaii Telescope (CFHT) in 
Hawaii and sporadically NARVAL on the 2 m Telescope 
Bernard Lyot (TBL) in southern France. ESPaDOnS and 
NARVAL are twin instruments able to collect stellar spectra 
spanning the whole spectral domain (from 370 to 1,000 nm) 
at a resolvin g power of 65.000, i n either circular or linear 
polarisation |Donati et alj l2~006cl . Donati et al., in prepa- 
ration). A total of 17 circular polarisation spectra were col- 
lected in 2006 February (9 ESPaDOnS spectra) and Novem- 
ber/December (7 ESPaDOnS and 1 NARVAL spectra); both 
data sets are referred to as Feb06 and Dec06 in the following. 



1 To ease comparison with BP Tau, we recall that V2129 Oph 
is 1.35 Mq star with a radius of 2.4 Rq, rotating with a pe- 
riod of 6.53 d and accreting mass at an esti mated rate o f abou t 
10~ s M Q yr _1 . According to the models of ISiess et"aT1 <200d) , 
V2129 Oph is no longer fully convective and hosts a small ra- 
diative core weighting about 0.1 Mq. Further information on 
V2129 Oph can be found in D07. 
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Table 1. Journal of observations. Columns 1—6 sequentially list the UT date, the instrument used, the heliocentric Julian date and UT 
time (both at mid-exposure), and the peak signal to noise ratio (per 2.6 kms" 1 velocity bin) of each observation (i.e., each sequence 
of four subexposures). Column 7 lists the rms noise level (relative to the unpolarized continuum level I c and per 1.8 kms -1 velocity 
bin) in the circular polarization profile produced by Least-Squares Deconvolution (LSD), while column 8 indicates the rotational cycle 
associated with each exposure (within each data set, and using the ephemeris given by Eq.[T]l. 



Date 


Ins t r u mc rrt 


HJD 


TTT 


+ 

L cxp 


^ /N 

O/ IN 


CT LSD 


Cycle 








(h.m.s) 


( s ) 




V-LU l c ) 




Feb 07 


ESPaDOnS / CFHT 


773.86074 


08:38:59 


4 X 600 


125 


5.4 


0.508 


Feb 08 


ESPaDOnS / CFHT 


774.72499 


05:23:37 


4 x 600 


130 


4.7 


0.622 


Feb 09 


ESPaDOnS / CFHT 


775.72732 


05:27:06 


4 x 400 


70 


11.0 


0.754 


Feb 10 


ESPaDOnS / CFHT 


776.72808 


05:28:20 


4 x 500 


120 


5.2 


0.885 


Feb 11 


ESPaDOnS / CFHT 


777.72195 


05:19:38 


4 x 600 


120 


5.5 


1.016 


Feb 12 


ESPaDOnS / CFHT 


778.72303 


05:21:20 


4 x 600 


120 


5.5 


1.148 


Feb 13 


ESPaDOnS / CFHT 


779.72259 


05:20:50 


4 x 600 


90 


8.5 


1.279 


Feb 14 


ESPaDOnS / CFHT 


780.72330 


05:21:59 


4 x 600 


70 


12.6 


1.411 


Feb 15 


ESPaDOnS / CFHT 


781.74836 


05:58:13 


4 x 600 


130 


5.4 


1.546 


Nov 29 


NARVAL/TBL 


1069.47610 


23:18:19 


4 x 900 


85 


7.3 


0.405 


Nov 30 


ESPaDOnS / CFHT 


1069.98455 


11:30:31 


4 x 900 


110 


3.5 


0.472 


Dec 05 


ESPaDOnS / CFHT 


1074.99539 


11:46:20 


4 x 900 


170 


2.7 


1.131 


Dec 06 


ESPaDOnS / CFHT 


1075.97733 


11:20:23 


4 x 900 


170 


2.8 


1.260 


Dec 07 


ESPaDOnS / CFHT 


1076.90044 


09:29:42 


4 x 900 


170 


2.9 


1.382 


Dec 08 


ESPaDOnS / CFHT 


1077.86802 


08:43:04 


4 x 800 


170 


2.8 


1.509 


Dec 09 


ESPaDOnS / CFHT 


1078.78979 


06:50:29 


4 x 600 


120 


1.2 


1.630 


Dec 10 


ESPaDOnS / CFHT 


1079.79637 


07:00:00 


4 x 600 


80 


6.6 


1.763 



Each polarisation spectrum consists of 4 individual subexpo- 
sures taken in different polarimeter configurations to remove 
all spurious polarisation signatures at the first order. 

Raw frames are processed with Libre ESpRIT, a fully 
automatic reduction package/pipeline developed and owned 
by JFD and installed both at CFHT and TBL. It auto- 
matically performs optimal extraction of ESPaDOnS and 
NARVAL unpolarised (Stokes 7) and circularly polarised 
(Stokes V) spectra following the procedure described in 
iDonati et alj (Il997l . Donati et al., in preparation). The ve- 
locity step corresponding to CCD pixels is about 2.6 kms -1 ; 
however, thanks to the fact that the spectrograph slit is 
tilted with respect to the CCD lines, spectra corresponding 
to different CCD columns across each order feature a differ- 
ent pixel sampling. Libre ESpRIT uses this opportunity to 
carry out optimal extraction of each spectrum on a sampling 
grid denser than the original CCD sampling, with a spectral 
velocity step set to about 0.7 CCD pixel (i.e. 1.8 kms - ). 
The peak signal-to-noise ratios (per 2.6 kms -1 velocity bin) 
achieved on the collected spectra (i.e., the sequence of 4 
subexposures) range between 70 and 170 depending on the 
weather and the exposure time. The full journal of observa- 
tions is presented in Table [T] 

Least-Squares Deconvolution (LSD: lDonati et al.lll997l ) 
was applied to all observations. The line list we employed 
for LSD is computed from an Atlas9 LTE model atmo- 
sphere (|Kurucj Il993l ) and corresponds to a K7IV spec- 
tral type (T e ff — 4,000 K and iogg = 3.5) appropriate 
for BP Tau. We selected only moderate to strong spectral 
lines whose synthetic profiles had line-to-continuum core de- 
pressions larger than 40% neglecting all non-thermal broad- 
ening mechanisms. We omitted the spectral regions within 
strong lines not formed mostly in the photosphere, such as 
the Balmer and He lines, and the Ca II H, K and infrared 
triplet (IRT) lines. Altogether, about 9,400 spectral features 
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Figure 1. LSD circularly-polarized and unpolarized profiles of 
BP Tau (top, bottom curves respectively) on 2006 Feb 11 (phase 
0.016). The mean polarization profile is expanded by a factor of 
25 and shifted upwards by 1.1 for display purposes. 



are used in this process, most of them from Fe I. Expressed 
in units of the unpolarised continuum level 7 C , the average 
noise levels of the resulting LSD signatures range from 2.7 
to 12.6xl0 -4 per 1.8 kms -1 velocity bin. 

All relevant information is summarised in Table[T] Rota- 
tional cycles E are computed from heliocentric Julian dates 
according to the ephemeris: 

HJD = 2453770.0 + 7.6E. (1) 
Note that the rotation cycle of BP Tau (7.6 d) was fully 
covered in Feb06, but only about 75% of it in Dec06. The er- 
ror on the period (0.1 d) translates into a phase uncertainty 
of 0.5 rotation cycle between the Feb06 and Dec06 data sets 
(separated by 298 d or 39 rotations). 
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Figure 2. Stokes / (panels 1 & 3) and Stokes V LSD profiles (panel 2 & 4) of photospheric lines of BP Tau (thick black line), for each 
observing night (top to bottom) of both Feb06 (left) and Dec06 (right) runs. The mean Stokes / LSD profile (averaged over the full set 
within each run, thin red line) is also shown to emphasise temporal variations. The rotational cycle associated with each observation is 
noted next to each profile; 3<r error bars are also shown to the left of each Stokes V profile. All profiles are plotted in the velocity rest 
frame of BP Tau. 



2.2 Zeeman signatures 

Zeeman signatures, featuring full amplitudes of about 0.5% 
of the unpolarised continuum level on average, are clearly 
detected in the LSD profiles of all spectra (e.g., see Fig. [TJ. 
The complete sets of LSD Stokes I and V profiles for both 
Feb06 and Dec06 runs are shown in Fig. [2] Note in particular 
that Zeeman signatures are significantly wider than unpo- 
larised profiles, suggesting that field strengths at the surface 
of BP Tau are rather strong. The corresponding longitudinal 
fields (com puted from the firs t order moment of the Stokes 
V profiles, iDonati et al.|[l997l ) are typically equal to a few 
hundred G; while they clearly vary with time, they keep the 
same (negative) sign throughout both runs (see Table [2]). 

We divided the LSD line list into 2 subsets, a blue 
and a red subset including lines bluer and redder than 
620 nm respectively; from those partial line lists, we derived 
LSD signatures (not shown) and corresponding longitudinal 
fields. We find that the longitudinal fields derived from the 
blue and red line subsets are slightly (about 15%) larger 
and smaller respectively than the average longitudinal fields 
derived from the original LSD signatures. This is similar 
(though much less extreme) than w hat was reported on the 
cool dwarf £ Boo (jPetit et al.ll2005l ); the origin of this effect 
is not clear. 

The longitudinal fields we de rive are significantly larger 
tha n the upper limits obtain ed bv lJohns-Krull et all |l999b) 
and IValenti fc Johns-Krull l|2004l ). This discrepancy is due 



to the different methods used to estimate longitudinal fields; 
the method used here (equivalent to the well-known "centre 
of gravity" technique) is known to be very robust for all fiel d 
configurations and strengths l|Landi Deel'Innocentil |2004|). 
while cross-correlation methods (used by I Johns-Krull et al.l 
1999b) are found to yield underestimates of both the longitu- 
dinal field and the corresponding error bar when Zeeman sig- 
natures are significantly broader than the unpolarised profile 
and feature a weak amplitude (as in the present case); long- 
term intrinsic variability of the large-scale magnetic topol- 
ogy may also partly explain the discrepancy. 

We also estimate the amount of veiling in the spectrum 
of BP Tau (at an average wavelength of 620 nm) by com- 
paring the equivalent widths of Stokes I LSD photospheric 
profiles between BP Tau and a standard star of similar spec- 
tral type (S Eri), using the same line list for both stars; the 
veiling parameter we derive varies between and 60% over 
the full span of our observations, i.e., co mparable to val- 
ues p reviously published in the literature (|Gullbring et alj 
Il998l ); the relative veiling variation over one single run is 
about ±15%, similar to previous reports by IValenti et al.l 
l|2003l ). 

Circular polarisation is also detected in most emis- 
sion lines, and in particular in the He I 587.562 nm (-D3) 
line and in the Ca 11 IRT lin es known as good trac- 
ers of magnetospheric accre tion l|johns-Krull et al1 Fl999b; 
IValenti fc Johns-Krulll I200II ) with medium magnetic sensi- 
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Table 2. Longitudinal magnetic field of BP Tau and associated error bars, as estimated from the LSD profiles (-BlsDi col. 2), the Ca II 
IRT emission core (Birt> col. 3), the He I 587.562 nm (D3) line (5jj c j D , col. 4), the He I 667.815 nm line (Bjj e j 667 , col. 5) and 
the 3 Fen lines of multiplet 32 around 500 nm (Bp e n , col. 6), Ho (B^ a , col. 7) and H/3 (Bjj^g, col. 8). Column 9 lists the veiling 
parameter r at each phase at an average wavelength of 620 nm (with typical error bars of about 0.01), defined as the relative difference 
in the equivalent width of Stokes I LSD photosphcric profiles between BP Tau and the standard star <5 Eri. Rotational cycles (col. 1) are 
computed according to Eq. [T] 



Cycle 




Birt 


B He 1 D 3 


^He I 667 


B Fe II 


%« 




r 




(kG) 


(kG) 


(kG) 


(kG) 


(kG) 


(kG) 


(kG) 




Feb06 run 


0.508 


-0.19 ±0.05 


0.53 ±0.07 


1.77 ±0.23 


5.32 ± 0.66 


2.84 ±1.10 


0.01 ±0.02 


0.61 ± 0.17 


0.22 


0.622 


-0.14 ±0.05 


0.44 ± 0.07 


1.44 ±0.18 


1.74 ±0.44 


1.64 ±0.61 


0.00 ±0.02 


0.43 ± 0.11 


0.29 


0.754 


-0.21 ±0.15 


0.64 ±0.12 


1.39 ±0.28 


3.93 ±0.50 


1.92 ±0.63 


0.02 ±0.03 


0.79 ± 0.34 


0.59 


0.885 


-0.33 ±0.07 


1.09 ± 0.07 


2.34 ±0.15 


4.36 ±0.33 


3.67 ±0.95 


0.08 ±0.02 


0.60 ±0.12 


0.35 


1.016 


-0.69 ±0.14 


1.57 ±0.08 


3.30 ±0.12 


5.58 ±0.25 


4.50 ±0.86 


0.23 ±0.02 


1.36 ± 0.15 


0.51 


1.148 


-0.63 ±0.13 


1.56 ± 0.09 


3.75 ±0.17 


8.17 ±0.47 


4.40 ±1.02 


0.12 ±0.02 


1.05 ± 0.14 


0.42 


1.279 


-0.50 ±0.16 


1.23 ± 0.12 


3.43 ±0.31 


7.88 ±0.88 


3.62 ±1.07 


0.14 ±0.03 


0.91 ± 0.25 


0.32 


1.411 


-0.19 ±0.15 


0.66 ± 0.14 


1.68 ±0.55 


7.45 ± 1.45 


3.93 ±1.67 


0.06 ± 0.03 


1.06 ± 0.41 


0.48 


1.546 


-0.31 ±0.07 


0.74 ±0.07 


1.85 ±0.19 


4.80 ±0.54 


2.67 ±0.91 


0.04 ±0.02 


0.36 ± 0.18 


0.21 


Dec06 run 


0.405 


-0.29 ±0.08 


0.75 ±0.08 


1.73 ±0.22 


3.96 ±0.43 


1.57 ±0.46 


0.06 ± 0.02 


0.15 ± 0.14 


0.39 


0.472 


-0.37 ±0.04 


0.92 ± 0.05 


2.36 ±0.12 


4.38 ±0.29 


2.64 ±0.58 


0.05 ±0.02 


0.80 ±0.09 


0.31 


1.131 


-0.28 ±0.03 


0.64 ± 0.04 


2.02 ±0.13 


4.95 ±0.35 


2.19 ±0.46 


0.02 ±0.02 


0.42 ± 0.07 


0.20 


1.260 


-0.53 ±0.05 


0.83 ± 0.04 


2.37 ±0.11 


3.71 ±0.22 


2.01 ±0.36 


0.04 ±0.02 


0.40 ±0.07 


0.39 


1.382 


-0.60 ±0.04 


0.72 ± 0.04 


1.85 ±0.09 


3.31 ±0.18 


2.64 ±0.49 


0.06 ±0.02 


0.52 ± 0.08 


0.41 


1.509 


-0.30 ±0.03 


0.93 ± 0.05 


2.62 ±0.12 


4.95 ±0.30 


3.10 ±0.81 


0.05 ±0.02 


0.93 ± 0.09 


0.20 


1.630 


-0.12 ±0.03 


1.17 ±0.07 


3.70 ±0.22 


7.78 ±0.72 


3.67 ± 1.00 


0.07 ±0.02 


0.46 ± 0.12 


0.06 


1.763 


-0.06 ±0.05 


0.83 ± 0.11 


2.28 ±0.38 


5.82 ± 1.34 


3.09 ±1.21 


0.07 ±0.03 


0.20 ±0.24 


0.00 



tivity (Lande factors of 1.0); clear signatures are also de- 
tected at all times in the He I 667.815 nm line (Lande fac- 
tor of 1.0) and in the Fe 11 492.393 nm, 501.842 nm and 
516.903 nm lines (multiplet 32, average Lande factor of 1.7). 
The complete sets of Stokes / and V profiles correspond- 
ing to the Ca 11 IRT0 and the He I D3 emission lines are 
shown in Figs. [3] and [4] respectively. The longitudinal fields 
(see Table [2} are significantly larger, and of opposite sign, 
than those derived from photospheric lines. They reach up to 
1.6 kG and 3.8 kG for the Ca 11 IRT and the He 1 D3 emission 
lines respectively, and up to 8.2 kG for the He I 667.815 nm 
line; longitudinal fields from the Fe 11 lines more or less re- 
peat (with larger error bars) those traced by the He I D3 line. 
Note that the field values we d erive from the He I D3 lin e 
are similar to those measured bv lJohns-Krull et al.l (ll999aT l. 
IValenti fe Johns-Krulll (|2004l ) and ISvmington et all |2005). 



3 ROTATIONAL MODULATION 

We now examine the temporal variations of the Stokes I and 
V profiles of photospheric and emission lines throughout the 
whole rotation cycle of BP Tau and demonstrate that these 
variations can be assigned to rotational modulation. 



2 Note that the 3 components of the Ca II IRT were averaged 
together in a single profile to increase S/N further. 



3.1 Photospheric lines and accretion proxies 

As for V2129 Oph (D07), the emission lines tracing accre- 
tion regions at the surface of the star (and in particular the 
Ca 11 IRT and the He I 587.562 nm and 667.815 nm lines) ex- 
hibit the strongest and simplest evolution over the timescale 
of our observations. The level of He I emission (e.g., Figs. [4] 
& O undergoes a progressive increase followed by a regu- 
lar decrease over th e rotation per i od, i.e ., very similar to 
what is reported by IValenti et all (|2003h . Moreover, com- 
paring, e.g., He 1 profiles at phases 0.508 and 1.546 on the 
Feb06 run demonstrates that the amount of emission more 
or less returns to its initial level once the star completed a 
complete rotation cycle (see Fig.0; the main part of the ob- 
served variability may thus be safely attributed to rotational 
modulation. 

Looking at the variability of the Ca 11 IRT emission pro- 
files in Feb06 brings the same conclusion (see Fig. [3}, even- 
though the amount of fluctuation - ±10% about the mean 
(i.e., from 18 to 22 kms -1 or 0.051 to 0.062 nm) - is smaller 
than that of the He I D3 line - almost a factor of 2 peak 
to peak (i.e., from 30 to 60 kms" 1 or 0.06 to 0.12 nm). 
The variation of the He I 667.815 nm line is even more ex- 
treme, reaching up to a factor of 3 peak to peak (from 6 
to 18 kms -1 or 0.013 to 0.040 nm); the Fe 11 emission lines 
vary in strength by about a factor of 2 (from 7 to 14 kms 
or 0.012 to 0.024 nm). Note however that the Call IRT 
emission profiles at phases 0.405 and 0.472 in the Dec06 run 
(both taken about 1 rotation cycle earlier than the bulk of 
the other Dec06 data) are significantly stronger than the 
average, and in particular stronger than what they should 
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Figure 3. The same as Fig. |2] but for the Ca II IRT lines. 
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Figure 4. The same as Fig.[2]but for the He I 587.562 nm (-D3) line. 
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have been in case of pure rotational modulation (given the 
observed profiles at phases 1.382 and 1.509). We attribute 
this to intrinsic variability, which we corrected by scal- 
ing down both profiles by 20% (Figs. [3] and [5] show pro- 
files/equivalent widths prior to this correction). Maximum 
emission is reached (in all lines simultaneously) at phases 
1.02 and 1.38 in the Feb06 and Dec06 runs respectively (see 
Fig. [5]). 

Stokes V profiles from emission lines strengthen the ev- 
idence that rotational modulation dominates the observed 
temporal variations. The Feb06 data are particularly clear 
in this respect; the shape of Zeeman signatures is indeed 
varying very smoothly with rotation phase (reaching max- 
imum strength around phase 1.02) and repeats well after 
one complete rotation cycle. The corresponding longitudi- 
nal field variations with rotational cycle, shown in Fig. [S] 
further confirm this impression. The Dec06 Zeeman data 
also show definite (though smaller) variability with rotation 
phase. The overall similitude between both data sets (with 
a smooth rise and decrease of the line emission and Zeeman 
signature amplitude) throughout the rotation period sug- 
gests that both line emission and Zeeman signatures were 
likely weak in the Dec06 phase gap during which no obser- 
vation were collected; we thus conclude that all the large 
magnetic regions have most likely been observed in Dec06 
despite the incomplete phase coverage. 

The longitudinal fields we measure, although grossly 
similar at both observing epochs, nevertheless exhibit a 
different rotational modulation in Feb06 and Dec06 (see 
Fig. [5} ; for instance, the longitudinal field curve shows only 
one maximum in Feb06, while it shows two distinct peaks 
in Dec06. It suggests that intrinsic variability likely dis- 
torted the large-scale field topology between both observ- 
ing epochs. Moreover, at both epochs, the longitudinal field 
curves are not symmetric with respect to the phase of maxi- 
mum line emission, suggesting that the field within accretion 
spots is not purely radial but rather slightly tilted with re- 
spect to the local meridional plane. In Feb06 for instance, 
the phases at which maximum longitudinal fields are ob- 
served are shifted with respect to the phase of maximum 
line emission (1.02) and differ for each line, with He I lines 
(and in particular the He I 667.815 nm line) showing max- 
imum phase shifts; it suggests that the tilt (with respect 
to the local meridional plane) of the magnetic field within 
accretion spots is larger in the formation region of He I lines. 

The Stokes V LSD photospheric profiles also show 
clear temporal variations (see Fig. [5} . Zeeman signatures in 
Feb06, tracing negative longitudinal fields throughout the 
whole rotation cycle, vary smoothly from small (e.g., phase 
0.622) to large (e.g., phase 1.016) amplitudes, then decrease 
again to their initial sizes and shapes (e.g., phase 1.546 and 
0.508). Similar behaviour is observed on Dec06 with maxium 
and minimum Zeeman signals around phase 1.38 and 1.76 
respectively. This variability strongly argues in favour of ro- 
tational modulation. Zeeman signatures are more complex 
than their emission line counterparts and trace a field of op- 
posite polarity; it suggests that emission and photospheric 
lines do not form over the same regions of the stellar surface. 



3 Other emission lines are mostly unaffected by this intrinsic vari- 
ability. 



A similar conclusion was reached in the case of V2129 Oph 
(D07). 

The Stokes / LSD profiles display only modest variabil- 
ity, apart from the changes in equivalent widths (±15% peak 
to peak within each run, up to 60% over both runs) that we 
attribute to veiling (see Table [2}. We find that statistically, 
veiling tends to be larger when emission lines are stronger; 
this correlation is howev er rather loose (as already pointed 
at bv lValenti et alj|2003l ). with strong veiling episodes some- 
times occuring when emission lines are weak (e.g., at phase 
1.411 in Feb06). Apart from veiling, variations in Stokes / 
photospheric profiles are also visible (though modest) and 
are reminiscent of bump-like signatures from dark spots 
travelling from the blue to the red side of the line profile 
(e.g., from phase 0.754 to phase 1.016 in Feb06, or from 
phase 1.260 to 1.382 in Dec06); at this point, and given the 
low v sin i of BP Tau, this assumption is no more than a 
speculation and surface imaging tools are needed to test it 
in more details. 

The average radial velocity of the photospheric lines is 
15.7 kms" 1 in Feb06 and 15.3 kms -1 in Dec06, suggesting 
an average radial velocity of about 15.5 kms -1 ; typical vari- 
ations of up to ±0.8 kms - are observed throughout each 
run. We assume in the following that this average position 
represents the heliocentric radial velocity of the stellar rest 
frame. The average width of photospheric lines (full width 
at half maximum of 15 kms -1 ) is mostly due to rotation 
(wsini = 9 kms -1 ). 

The Ca II IRT emission core is centred at 16.5 kms - 
on average, i.e., redshifted by about 1 kms -1 with respect 
to the photospheric lines. As for photospheric lines, varia- 
tions of up to ±0.8 kms -1 are observed throughout each 
run, with a smooth dependance with rotational phase (see 
Fig. [5}. The width of the central emission core (full width 
at half maximum of 21 kms -1 ) is comparable to the rota- 
tional broadening of the star. The Fe II emission lines are 
centred at 17.3 kms -1 on average and feature a full width 
at half maximum of about 24 kms -1 . Velocity variations are 
stronger, reaching up to ±2 kms -1 throughout both runs. 

Both He I emission lines exhibit higher redshifts rela- 
tive to the stellar rest frame than the Ca II IRT emission 
core and lie at 23 kms -1 on average, i.e., about 7.5 kms -1 
redward of the photospheric lines. Larger variations of the 
line centroids (of ±2 — 3 kms -1 ) are also observed through- 
out each run; the phase dependance is rather smooth and 
repeats well between both lines, suggesting that it is caused 
by rotational modulation (see Fig. [SJ ■ Since maximum He I 
emission occurs roughly half-way through the blue to red 
migration of the He I emission lines (i.e., at phase 1.00 and 
1.40 for the Feb06 and Dec06 runs respectively), we specu- 
late that this is likely where the parent accretion spots are 
located. While the He I 667.815 nm line is roughly as wide 
as the Ca II emission cores (full width at half maximum of 
22 kms -1 ), the He I D3 line is twice as wide (full width at 
half maximum of 44 kms -1 ) as a result of being a composite 
profile of 6 different transitions. 

A broad Ca 11 emission component (with a full-width at 
half ma ximum of about 200 k ms -1 , similar to what shown 
by, e.g.. lGullbring et al.l li996 in their Fig. 6) is also present 
(though not shown here) in some of our spectra, mostly 
those from Dec06. Its red wing more or less anticorrelates 
with the narrow central component, featuring an emission 
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Figure 5. Equivalent widths (top), longitudinal fields (with ±1<t error bars, centre) and line centroids (bottom, the stellar velocity rest 
frame being at 15.5 kms ) as derived from the Call IRT (red circles), the He I D3 (green squares) and the He I 667.815 nm (blue 
triangles) emission lines as a function of rotational cycle, for both Feb06 (left) and Dec06 (right) runs. Relative error bars on equivalent 
widths (top panels) are smaller than the symbol size. Note the logarithmic vertical scale for the top and centre panels. 



bump around phases of minimum core emission (phase 1.63 
in Dec06) and a relative dip (with respect to the mean pro- 
file) around phases of maximum core emission (phase 1.26- 
1.38 in Dec06). As it correlates well with the red wing of 
Balmer lines, we attribute it to disc material free-falling 
along accretion funnels (see below). This broad component 
is apparently subject to a higher level of intrinsic variabil- 
ity than the narrow component. Stokes V signatures being 
entirely confined to the central narrow emission component, 



we do not consider the broad component any further in the 
present study and simply removed it from all Ca n profiles 
(and in particular from those shown in Fig. [3J. No similar 
broad emission component is visible in He I and Fe 11 lines. 

The widths and redshifts of the Ca II, Fe 11 and He I Zee- 
man signatures are similar to those of the unpolarised emis- 
sion profiles. However, while the Stokes V signatures from 
the Ca 11 and Fe 11 emission lines are roughly antisymmetric 
(with respect to the line centre), those from He I lines show 
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significant departures from antisymmetry, with the biue iobe 
being both narrower and deeper than the red lobe. It sug- 
gests that Zeeman signatures from He I lines form in regions 
featuring significant velocity gradients, whereas Zeeman sig- 
natures from Ca II and Fe II lines form in regions that are 
almost at rest. This is in good agreement with the observed 
line velocity redshifts, much larger for the He I lines than for 
the Ca II and Fe II lines. We also note that both He I lines 
trace larger toroidal fields than the Ca II and Fe II lines (see 
above) . 

Most observations reported here are similar to what 
V2129 Oph exhibits (D07), suggesting that the same mod- 
elling strategy can be employed. We thus assume, as for 
V2129 Oph, that emission lines comprise two physically dis- 
tinct component^. We attribute the first of these, the accre- 
tion component, to localised accretion spots at the surface 
of the star whose visibility varies as the star rotates, giving 
rise to rotational modulation of the emission (both in in- 
tensity and radial velocity). A second, chromospheric com- 
ponent is more or less evenly distributed over the surface 
of the star, producing a time-independent emission compo- 
nent. We further assume that the chromospheric component 
is mostly unpolarised while Stokes V signatures arise only 
in the accretion component. With the relative strengths of 
both components varying from line to line, we expect dif- 
ferent lines to yield different longitudinal fields; moreover, 
we expect lines showing higher levels of variability to yield 
stronger longitudinal fields and larger velocity variations, in 
agreement with what we observe on BP Tau. This model is 
further detailed in Sec. |4] 

3.2 Balmer lines 

Ha and H/3 lines exhibit strong emission with average equiv- 
alent widths of 4,850 and 1,210 kms -1 , (i.e., 10.6 and 
2.0 nm) respectively. Maximum emission is reached at phase 
0.89 and 1.38 in the Feb06 and Dec06 runs respectively, i.e., 
at roughly the same time as all other emission lines (see 
above). It suggests that rotational modulation dominates 
the observed varibility, reaching about ±20 — 25% for H/3 
and ±10% for Ha. The variations are however not as smooth 
as those of the accretion proxies described above (e.g., be- 
tween phase 0.62 and 0.88 in Dec06), suggesting tha Balmer 
lines are subject to a higher level of instrinsic variability. 

Both Ha and H/3 exhibit Stokes V signatures (not 
shown here) similar in shape to those of V2129 Oph (D07), 
with significant departures from antisymmetry. The corre- 
sponding longitudinal fields (listed in Table [2j|, are peaking 
at about 0.23 and 1.4 kG for Ha and H/3 respectively; they 
are smaller than, and correlate reasonably well with, those of 
other emission lines, suggesting that rotational modulation 
dominates the observed variations. The weaker longitudinal 
fields indicate that the circularly polarised signal from the 
accretion spots is more diluted (with the unpolarised light 
from the chromosphere and wind, and with weakly polarised 
contributions from the lower-field outer regions of accretion 

4 These two modelling components should not be confused with 
the narrow and broad components of the Ca II emission lines men- 
tionned above, and are both used to describe the narrow line 
emission core only. 



funnels) in Balmer lines than in the other accretion proxies; 
we speculate that most of the intrinsic variability observed in 
Ha and H/3 is due to this unpolarised chromospheric/wind 
component or to that from the outer regions of accretion 
funnels. 

We also note that the unpolarised profile of Ha and 
H/3 are much broader (with full widths at half maximum of 
about 240 kms -1 ) than their Zeeman signatures (70 kms -1 
wide). Moreover, both unpolarised lines are significantly 
blueshifted with respect to the stellar velocity rest frame 
(by as much as 15 kms" 1 ) while their Stokes V signatures 
are slightly redshifted (by a few kms -1 , as Ca II and Fe II 
emission lines). It confirms that Balmer lines contain sev- 
eral emission components, with at least a slightly redshifted 
one forming in accretion spots and producing the observed 
Stokes V signatures, and a second (dominant) blueshifted 
one due to a chromosphere and/or a wind. This is all very 
reminiscent of what is observed on V2129 Oph (D07). 

Balmer lines collected in Dec06 also include a conspic- 
uous high-velocity component in their red wings (between 
±100 and ±300 kms" 1 , and mostly visible in H/3). This 
component varies roughly in phase opposition with the rest 
of the line and shows up as a hump around phase 1.63 (i.e., 
when the overall line emission is small) and as a relative 
dip (with respect to the average profile) around phase 1.26- 
1.38 (i.e., when line emission is large). This is readily vis- 
ible on the corresponding H/3 autocorrelation matrix (see 
Fig. 0. By analogy with V2129 Oph (D07), we attribute 
it to disc material free falling onto the star within accre- 
tion funnels. Given the limited agreement between profiles 
taken almost one complete rotation cycle apart (e.g., phase 
0.405 and 1.382, or phases 0.472 and 1.509), we conclude 
that component is apparently subject to a high level of in- 
trinsic variability. This component is also present, though 
much fainter, in the Feb06 data set; it correlates well with 
the red wing of the broad Ca II IRT emission component 
discussed above, thus suggesting a common origin. 



4 THE ACCRETION SPOTS AND MAGNETIC 
TOPOLOGY 

Now that the rotational modulation of Stokes / and V pro- 
files of photospheric and emission lines is firmly established, 
we aim at modelling the observed modulation and derive 
from it maps of the accretion spots and the magnetic field 
at the surface of BP Tau. As in D07, we only use the LSD 
photospheric and Ca II line profiles in this study. Modelling 
He I lines (and in particular the non-antisymmetric Stokes V 
profiles, see Fig. 3]) is rather uncertain; it indeed requires ad- 
ditional independent information on the velocity fields and 
gradients within the line formation regions that we usually 
do not have and cannot easily access. 

4.1 Model description 

The model we use here is directly inspired from that of D07. 
It assumes that Stokes V photospheric and Ca II emission 
lines form in different regions of the stellar surface, reflecting 
the fact that accretion spots (where emission lines mostly 
form) coincide with dark cool spots at photospheric level 
(to which photospheric lines are mostly insensitive). This 
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Figure 6. Temporal variations of the Ha (panels 1 & 3) and H/3 (panels 2 & 4) profiles of BP Tau (thick black line), for each observing 
night (top to bottom) of both Feb06 (left) and Dec06 (right) runs. The mean profile (averaged over the full set within each run, thin red 
line) is also shown to emphasise temporal variations. The rotational cycle associated with each observation is noted next to each profile. 
All profiles are plotted in the velocity rest frame of BP Tau. 
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Figure 7. H/3 autocorrelation matrix computed from the Dec06 
data. White and dark brown indicate perfect correlation and an- 
ticorrelation respectively. 



is compatible with reports that BP Tau looks darkest in 
optical photometry when accretion spots are most visible 
i|Valenti et al.ll2003l ). 

We assume that accretion lines can be described with 
the two component model introduced in Sec. [3] combin- 
ing a chromospheric emission component evenly distributed 
over the star with an additional emission component con- 
centrated in local accretion spots. The emission-line Zee- 
man signatures are assumed to be associated with the accre- 
tion component only. In this model, rotational modulation 
of emission lines results from the accretion regions being 
carried in and out of the observer's view as the star rotates, 
the excess emission received from these regions being pro- 
portional to their projected area (as observations suggest). 

We introduce a vector magnetic field B and a local ac- 
cretion filling-factor /, describing the local sensitivity to ac- 
cretion and photospheric line^f). For / = 0, the local area 
on the protostar's surface contributes fully to photospheric 
lines and generates no Ca n excess emission (and only un- 
polarised chromospheric Call emission); for / = 1, the lo- 
cal area does not contribute at all to photospheric lines and 
produces the maximum amount of excess Ca II emission and 
polarisation (in addition to the unpolarised chromospheric 



model was fairly successful at reproducing observations of 
V2129 Oph. Given the overall similarity of the present data 
sets with that of V2129 Oph (see above), we speculate that 
it should be adequate for BP Tau as well; in particular, it 



5 Note that the local accretion filling-factor we define here is 
different from the usual accretion filling factor of the cTTS liter- 
ature, i.e., the relative area of the total stellar surface covered by 
accretion spots. 
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Ca II emission). Spectral contributions for intermediate val- 
ues of / are derived through linear combinations between 
the f — and / = 1 cases. 

We introduce one difference with the model of D07. 
We assume that, in each local cell of the stellar surface, 
only a fraction ip of the cell area contains magnetic fields 
(whose local strength thus equals B/ip). We further assume 
that, within each surface cell, the accretion filling factor de- 
fined above only applies to the magnetic portion of the cell, 
while the non-magnetic portion corresponds to / = 0. We 
finally assume (for simplicity) that ip is constant over the 
whole stellar surface. In this context, the model of D07 cor- 
responds to the specific case where ip = 1. Introducing ip 
brings the possibility of fitting magnetic fluxes (B) and mag- 
netic strengths (B/ip) independently, both having different 
effects in spectral lines of stars with low to moderate v sin i. 
While the magnetic strength is mainly derived from fitting 
the far wings of the Zeeman signatures, the magnetic flux is 
obtained essentially by adjusting the amplitude of the Stokes 
V profile^. We consider ip, called the magnetic filling fac- 
tor, as a free parameter that we optimise by minimising the 
amount of information in the magnetic image (see below) . 

In this context, the local synthetic photospheric and 
emission Stokes 7 and V lines profiles emerging from each 
grid cell (noted respectively 7 P , V p , 7 C and V c ) are given by: 



Ic 

Ve 



(l - /)v»/ m + (l - 
(i - f)i>v m 

Ic + Mh - 1) 



(2) 
(3) 
(4) 
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where 7 m and 7 q are the local Stokes I photospheric pro- 
files corresponding to the magnetic and non-magnetic areas, 
Vm the local Stokes V photospheric profile corresponding 
to the magnetic areas, 7 C and 7 a the local Stokes I profiles 
corresponding to the chromospheric and accretion emission 
components, and V a the local Stokes V profile corresponding 
to the accretion emission component. 

We obtain both B and / by fitting the correspond- 
ing synthetic Stokes V and 7 spectra to the observed Zee- 
man signatures and unpolarised LSD profiles of photospheric 
lines and Ca II emission cores. The code we use for fitting B 
and / is adapted from the stellar surface magne tic imaging 
code of lDonatil (|200lh and IDonati et alj (|2006rJ ). The field 
is decomposed in its poloidal and toroidal components (with 
the poloidal field split itself between its radial and non ra- 
dial contributions) and described as a spherical harmonics 
expansion, whose coefficients (ae, m , (3e, m an d "ie.m for the 
radial, non-radial poloidal and toroidal components respec- 
tiv ely, where t and m de note the modes orders and degrees, 
see IDonati et al. 2006b for more details) are unknowns in 
the fitting procedure. The remaining free parameters are 
the values fj of / over the surface of the star (divided into 
a grid of thousands of small surface pixels indexed with j). 

The inversion problem being ill-posed, we use the princi- 
ples of maximum entropy image reconstruction to make the 



6 Magnetic distortions on Stokes I profiles remain small, at least 
on the average photospheric line considered here, much smaller 
in particular than the observed rotational modulation described 
in Sec. [3] 



solution unique (|Skilling fc Brvanll 19841 ). The entropy func- 
tion, computed from fj, ae, m , Pe.m and 7£, m , allows us to 
select the image with minimum information (maximum en- 
tropy) given a predetermined quality of the fit to the data, 
usually set to a reduced chi-square x? of about 1. By fit- 
ting only at,m an d Pe,m, we have the possibility of trying to 
force the solution towards a purely poloidal field; similarly, 
by weighting odd or even coefficients very heavily in their 
contribution to the entropy function, we also have the option 
of driving the recovered field topology towards symmetry or 
antisymmetry with respect to the centre of the star. An- 
other significant advantage of this new imaging method is 
that it suffers much smaller crosstalk between field compo- 
nents than the original one IjDonati fc Brownlll997l ; lDonatil 
l200lfi . 

We use Un n o-Rac hkovsky's equations (e.g., 
Landi Pegl'Innocentl! |2004) to model the LSD profiles 
of local photospheric lines (i.e., 7 q , I m and V m ), whose 
equivalent wavelength and Lande factors are set to 620 nm 
and 1.2 respectively. We adjust the free Milne-Eddington 
model parameters by fitting Unno-Rachkovsky's equations 
to an unpolarised LSD profile of the slowly rotating stan- 
dard star S Eri whose spectral type is very similar to that 
of BP Tau. We then derive both the rotational broadening 
v sin i and the radial velocity u ra d of BP Tau by fitting our 
series of LSD Stokes 7 profiles and selecting the values that 
minimise the image information content (for a given quality 
of the fit to the data). We find usini = 9.0 ±0.5 kms -1 and 
«rad = 15.5 ± .5 kms" , in good agreem ent with previous 
estimates fe.g.. Ijohns-Krull et al.|[l999bl ). 

The local Ca II Stokes I chromospheric and accretion 
profiles I c and J a are both described (in the absence of mag- 
netic fields) with a simple gaussian centred at 850 nm, whose 
full width at half maximum (about 15 kms -1 ) is derived 
from a fit to the observed Ca II emission profiles of BP Tau 
(assuming the rotational broadening obtained in the previ- 
ous modelling step). The effect of magnetic fields on 7 a and 
the corresponding Stokes 7 and V profiles are modelled using 
Unno-Rachkovsky's equations and assuming a unit Lande 
factor. The equivalent width ratio of the accretion emission 
component 7 a to that of the quiet chromospheric compo- 
nent 7 C is set to A/ip (the division by ip ensuring that the 
total amount of local emission, and hence the derived ac- 
cretion filling factor map, does not depend on ip); similar- 
results being obtained for other values, e.g., 3/ip or 5/tp. 
The equivalent width eb of the quiet chromospheric emis- 
sion component 7 C is a free parameter that we optimise by 
minimising the information content of the magnetic image 
(see below); increasing eb will force the imaging code to de- 
crease the relative fraction of accretion emission with respect 
to that of chromospheric emission (and vice versa), so that 
the resulting synthetic emission profiles match the observed 
ones. 

4.2 First attempt 

The first attempt consists in fitting, for both data sets, LSD 
Stokes 7 photospheric profiles only. In addition to provide 
an accurate estimate of both usini and ?; ra d (see above), 
this step allows us to check the validity of our assumption 
that the shape of photospheric lines is mostly distorted by 
dark spots at the surface of the star and that the observed 
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Figure 8. Stokes I LSD profiles of photospheric lines of BP Tau 
(thick black line) along with the maximum entropy fit (thin red 
line) to the data, for both Feb06 (left) and Dec06 (right) runs. 
The rotational cycle of each observation is noted next to each 
profile. 



profile variability is compatible with rotational modulation. 
It also enables us to check a second modelling assumption, 
i.e., that the putative dark photospheric spots at the surface 
of the star are indeed more or less coincident with the accre- 
tion regions identified from the intensity and radial velocity 
modulation of emission lines (see Sec. [3}. 

We start by removing veiling from all LSD Stokes I pho- 
tospheric profiles, i.e., by scaling them to the same equiv- 
alent width; modelling veiling as part of the imaging pro- 
cess appears undesirable at this stage given that veiling 
only weakly correlates with the other parameters we aim at 
modelling (and in particular the Zeeman signatures and the 
shape and strength of photospheric and emission line pro- 
files). The fit we obtain (see Fig. [8} matches the data at a 
S/N level of about 500 and yields a Xr improvement of about 
a factor of 2 with respect to synthetic profiles corresponding 
to an unspotted star. Most of the main profile distortions, 
like for instance the profile asymmetries at phase 1.260 and 
1.382 in the Dec06 data set (see Sec. [3] and Fig. [2}, are 
reproduced by our model, leading us to conclude that our 
series of Stokes I profiles of BP Tau are compatible with 
rotational modulation induced by dark surface spots. The 
reconstructed images (see Fig. show one main spot at 
each epoch, covering in both cases about 2% of the total 
stellar surface. 

These images should not be viewed as true brightness 
images of the photosphere of BP Tau; given the rather low 
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Figure 9. Maps of the local surface brightness (relative to that 
of the quiet photosphere) on BP Tau, for both Feb06 (top) and 
Dec06 (bottom) runs. The star is shown in flattened polar pro- 
jection down to latitudes of —30°, with the equator depicted as 
a bold circle and parallels as dashed circles. Radial ticks around 
each plot indicate the phases of observation. 

v sin i (by Doppler imaging standards) and the subsequently 
limited spatial resolution (about 0.1 cycle at the rotation 
equator), only the largest non-axisymmetric spots are re- 
constructed here. In particular, polar spots such as that de- 
tected on V2129 Oph (producing no rotational modulation) 
are hardy detectable on BP Tau through Stokes I profiles 
only. The spots we recover should thus be seen as the most 
reliable features without which the distortions and variabil- 
ity of unpolarised profiles cannot be interpreted in terms of 
rotational modulation. 

We find that the reconstructed dark spots are reason- 
ably close to where we expect accretion regions to concen- 
trate, i.e., at high latitudes - given the low amplitude of ve- 
locity variations of emission lines (see Fig. [5| - and at phase 
1.02 (in Feb06) and 1.38 (in Dec06) - at which emission 
lines show maximum emission and reach a median position 
in their blue to red transit (see Figs. [3j [4] and [5} ■ They ap- 
parently lag slightly behind (by about 0.05 to 0.10 rotation 
cycle) the accretion spots themselves. Although the spatial 
resolution is likely too low to ascertain the latter point, it is 
nevertheless enough to confirm our basic working assump- 
tions. 

At this stage, one may argue that that the rotational 
broadening we derive in this process is overestimated as it 
neglects line broadening induced by potentially strong mag- 
netic fields at the surface of BP Tau; moreover, some of the 
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Figure 10. Stokes V LSD profiles of photosphcric lines of BP Tau 
(thick black line) along with the maximum entropy fit (thin red 
line) to the data, for both Feb06 (left) and Dec06 (right) runs. 
The rotational cycle of each observation and 3cr error bars are 
also shown next to each profile. 



profile distortions we trace and attribute to dark spots could 
be due to magnetic fields instead. We however think that is 
not the case, for at least two main reasons. Firstly, the v sin i 
we deri ve is fully compati b le (and even slightly lower) than 
that of IJohns-Krull et al.l (Il999bh . demonstrating that our 
estimate is obviously not an upper limit. Secondly, we would 
expect to see significant changes with rotational phase in 
the width of the LSD unpolarised profile (as we do for in- 
stance in_chenrically p eculia r stars hosting strong magnetic 
fields. lLandstreetl e.g.. [l988l ) since the magnetic topology of 
BP Tau is apparently not fully axisymmetric; no such effect 
is detected. 

Yet, the magnetic fields detected on BP Tau, reaching 
strengths of up to 8 kG (see Table [2}, are in principle strong 
enough to produce straightforwardly visible broadening of 
unpolarised lines. So why do we see no such signatures in 
our series of Stokes I LSD profiles? This is because very in- 
tense fields actually concentrate where accretion is strongest 
and thus hide in very dark photospheric regions; for this rea- 
son, they mainly show up in emission lines and almost disap- 
pear from optical photospheric lines, whose widths therefore 
remain mostly unaffected. This interpretation (mostly spec- 
ulative at this point) is confirmed below. 
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Figure 11. The same as Fig.|8]but for the Ca II emission profiles 
of BP Tau. 
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Figure 12. The same as Fig. llOl but for the Ca II emission profiles 
of BP Tau. 
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Figure 13. Maps of the local accretion filling factor / on BP Tau, 
for both Feb06 (top) and Dec06 (bottom) runs. Note that these 
maps also depict the reconstructed photospheric brightness (equal 
to 1 — ipf), our model being based on the assumption that dark 
photospheric spots and accretion regions coincide. These maps, 
obtained when modelling all data simultaneously, can be com- 
pared to those of Fig. Q|] derived from fitting Stokes / profiles of 
photospheric lines only. 



4.3 Complete modelling 

In this second step, we carry out the complete modelling 
by fitting, for each observing epoch, a magnetic topology 
and an accretion filling-factor map to all Stokes I and V 
profiles of photospheric and emission lines simultaneously. 
The remaining free parameters to adjust in this process are 
the magnetic filling factor ip (arbitrarily set to 1 in the first 
modelling step) and the equivalent width of the quiet 
chromospheric emission component. We find that ip needs 
to be significantly smaller than 1 for LSD Stokes V profiles 
of photospheric lines to be fitted down to noise level, in 
particular in the far line wings. 

Given the available spatial resolution in the data, spher- 
ical harmonics expansions describing the magnetic field were 
truncated at £ — 10. Moreover, the field reconstruction was 
oriented towards antisymmetric (rather than symmetric) 
magnetic topologies (about the center of star), by favour- 
ing spherical harmonics terms with odd I values (see D07); 
dominantly antisymmetric field configurations are indeed 
the only ones capable of yielding (through their dipolar 
component in particular) the exclusive high-latitude anchor- 
ing of accretion fun nels that we observe for BP Tau (e.g., 
iGreeorv et ai1l2006l ). 



The value that minimises the amount of energy in the 
reconstructed field given a unit x? fit to the data is ip = 0.25. 
We also find that eb should be about 60% of the average 
observed Ca II line emission, the remaining 40% being pro- 
duced by accretion spots and modulated by rotation. The 
final fits to the data, corresponding to a unit x?, are shown 
in Fig. [10] for the Stokes V LSD photospheric profiles, in 
Fig.[TTJfor the Stokes / Ca II emission lines and in Fig. ll2l for 
the Stokes V Ca II emission lines; the new fit to the Stokes / 
LSD photospheric profiles (not shown) is only slightly worse 
than that of Fig. [8] The reconstructed accretion filling factor 
and magnetic flux maps are shown in Figs. ll3l and ll4l respec- 
tively. Since accretions regions are assumed to concide with 
dark photospheric spots in our model, Fig [13] also provides 
a description of the reconstructed photospheric brightness. 

The updated v sin i estimate we derive (as a by-product) 
from this new modelling is virtually identical to that ob- 
tained in the previous modelling step (involving no magnetic 
fields) ; it confirms in particular that magnetic distortions on 
Stokes / profiles are weak (at least on the average line that 
we consider here). Attempts at fitting the data (and in par- 
ticular the Stokes V sets) assuming different values of the 
rotation period confirm that optimal results are obtained for 
a rotation period close to the nominal value of 7.6 d (within 
about 0.5 d). 

The reconstructed accretion regions (and dark photo- 
spheric spots) are mostly located over the polar regions of 
BP Tau; each spread out on about 8% of the star, actually 
covering only 2% of the total surface given than only one 
fourth of each surface pixel is subject to accretion and hosts 
magnetic fields {ip = 0.25, see model description above). 
While not identical to the initial brightness maps of Fig. [9] 
(derived from fitting Stokes / profiles of photospheric lines 
alone in our first modelling attempt), the new accretion fill- 
ing factor (and photospheric brightness) maps that we now 
obtain from fitting the full data set are nevertheless similar; 
in particular both feature mid-latitude appendages located 
at more or less the same phase as the brightness spots of 
Fig. [5] We suspect that some of the discrepancy between 
both sets of maps is due to the fact that the assumptions 
underlying our model (e.g., dark spots overlapping accre- 
tion regions) are too simple; however, since they both yield 
acceptable fits to the Stokes / LSD profiles, we conclude 
that their differences likely reflect the limited spatial reso- 
lution available in the spectra of BP Tau. Our result con- 
firms that the magnetic broadening of Stokes I LSD pro- 
files is small, thereby confirming the preliminary conclusions 
reached above. 

The magnetic topologies we recover at both epochs are 
similar, apart from a phase shift of about 0.25 rotation cycle. 
They both include one main positive radial field feature close 
to the magnetic pole, where the magnetic flux reaches 3 kG 
and the field strength up to 12 kG (since ip — 0.25). The av- 
erage field strength over the accretion regions is about 9 kG, 
in agreement with the highest values of longitudinal fields 
traced by emission lines. In non accreting regions, the aver- 
age magnetic flux is about 1.2 kG; the field is dominantly 
meridional and pointing away from the observer, matching 
the constantly negative longitudinal fields traced by photo- 
spheric lines. 

We find that the magnetic topologies we recover are 
dominantly poloidal, with only about 10% of the magnetic 
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Figure 14. Magnetic topologies of BP Tau in Feb06 (top) and Dec06 (bottom), reconstructed from a simultaneous fit to the complete 
series of Stokes I and V profiles of photospheric lines and Ca II emission cores. The three components of the field in spherical coordinates 
are displayed (from left to right), with magnetic fluxes labelled in G. 



energy concentrating in the toroidal component. In particu- 
lar, the toroidal field of BP Tau is definitely smaller (in terms 
of fractional magnetic energy) than that of V2129 Oph, 
which gathered about 20% of the reconstructed magnetic 
energy (D07). This toroidal component seems nonetheless 
real; while fitting the Stokes V data at unit Xt with a purely 
poloidal field is possible, the magnetic field we recover is sig- 
nificantly stronger suggesting that a purely poloidal solution 
is far less likely (according to maximum entropy principles). 

The spherical harmonics terms dominating the recov- 
ered magnetic topology corresponds to a slightly tilted 
dipole, with about 50% of the magnetic energy concentrat- 
ing in I = 1 modes; the corresponding dipole strength is 
about 1.2 kG. The second dominant term is a slightly tilted 
octupole, with about 30% of the magnetic energy gather- 
ing in £ — 3 modes; the corresponding octupole strength 
is about 1.6 kG0. This makes BP Tau fairly different from 
V2129 Oph, whose field comprises a dominantly octupolar 
large-scale magnetic field and only a small dipole compo- 
nent (D07). All remaining spherical harmonics terms of the 
large-scale (i.e., £ < 10) poloidal field of BP Tau contain 
altogether no more than 10% of the total magnetic energy. 
Note that our data does not exclude the potential presence 



7 While the dipole and the octupole are both tilted by about 10° 
with respect to the rotation axis, they are not tilted towards the 
same phase so that their poles do not coincide. 



of small-scale magnetic features (producing no detectable 
polarisation signatures) at the surface of BP Tau. 

Although small, the intrinsic differences between the 
two reconstructed magnetic topologies of BP Tau are ap- 
parently real and can be traced directly to differences in 
the data sets themselves; Zeeman signatures in photospheric 
lines around phases of maximum longitudinal fields (i.e., 
at cycles 1.016 and 1.148 in Feb06 and 1.260 and 1.382 in 
Dec06) are indeed different in shape and undergo different 
temporal evolution. These changes are however fairly limited 
and do not affect the gross characteristics of the large-scale 
field. 



5 MAGNETOSPHERIC ACCRETION AND 
CORONA 

As for V2129 Oph, we propose here an illustration of how 
accretion may proceed between the inner disc and the sur- 
face of the star. To do this, we extrapolate the reconstructed 
magnetic field over the whole magnetosphere, assuming the 
3D field topology is mainly potential and becomes radial be- 
yond a certain magnetospheric radius ta from the star (to 
mimic the openin g of the largest magnetic loops under the 
coron al pressure, [jardine et all 12003. [20061 : iGregorv et al.l 
2006). In non-accreting stars, this distance is usually as- 
sumed to be smaller than or equal to the corotation radius 
(rc) at which the Keplerian orbital period equals the stellar 
rotation period. In cTTSs however, the magnetic field of the 
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protostar is presumably clearing out the central part of the 
accretion disc and extending as far as the inner disc rim. 

This is of course only an approximation. In particular, 
the strong plasma flows linking the disc to the stellar surface 
as a result of mass-accretion will likely prevent the field from 
being potential by building up a strong azimuthal compo- 
nent in the magnetosphere. Magnetic tracers at the base of 
accretion funnels (and in particular He I emission lines, see 
Sec. [3) suggest that this is indeed the case. We nevertheless 
use this approach as a first step; further detailed simulations 
are postponed to forthcoming papers. The magnetospheric 
maps we derive for BP Tau (see Fig. [15] for two possible 
values of ta at epoch Feb06) show that the field topology is 
complex close to the stellar surface, but dominated by more 
extended open and closed field lines at larger distances. 

From these maps, we can estimate where the accretion 
funnels are located and where they are anchored at the sur- 
face of the star. We do this by identifying those magneto- 
spheric field lines that are able to accrete material from the 
disc, i.e., that link the star to the disc and intersect the ro- 
tational equator with effective gravity pointing inwards in 
the c o-rotating frame of reference (e.g. iGregorv et al.ll2006l . 
120071 ). We find that accretion spots concentrate at high lat- 
itudes if we assume ta = 7.5 i?* (see Fig. [15) . For values 
of 7" a smaller than about 3.5 i?*, equatorial accretion spots 
start to form; most of the accretion occurs onto the equator 
when r\ < 2.5 R+. 

Given the results of Sec. [4] demonstrating that accre- 
tion spots are located at high latitude, we conclude that va 
is at least equal to about 4 R* and may extend as far out 
as re, equal to about 0.07 AU or 7.5 i?* for BP Tavfl Our 
modelling also demonstrates that accretion spots concen- 
trate at phase of maximum line emission (0.0 in Feb06 and 
0.4 in Dec06, see Fig. [13); this is further confirmed by the 
transient absorption appearing in the red wing of Balmer 
lines at these phases (see Sec. [3). When assuming that the 
accretion disc lies within the equatorial rotation plane, the 
northern accretion hot spot that our extrapolated magne- 
tosphere model predicts (located at longitude 90° or phase 
0.75 in Feb06, see Fig. [I6j left panel) does not match the 
observations. Assuming that the accretion disc is magneti- 
cally warped around ta and locally lies within the magnetic 
equator of the large-scale dipole component produces a bet- 
ter agreement with observations (see Fig. 1161 right panel). 



6 DISCUSSION 

The new spectropolarimetric data we collected on BP Tau 
at 2 different epochs enabled us to obtain a realistic model 
of the large-scale magnetic field on this prototypical cTTS. 
Comparing our results with those recently obtained for a 
more massive cTTS (V2129 Oph, D07) suggests new hints 
on how magnetic fields are produced in young low-mass stars 



8 Assuming that the s poradic fluctuation s in the period of photo- 
metric variations (e.g. . S imon et al Jl990h is due to contamination 
by optical light from the inner disc rim (presumably modulated 
on a timescale corresponding to the local Keplerian period), it 
would imply that the radius of the inner disc, equal to Ta, also 
fluctuates with time, from about 6.5 to 8.5 i?* typically. 



and how they make young stars interact with their accretion 
discs. 



6.1 Dark spots and accretion regions 

As for V2129 Oph, we observe that photospheric lines and 
the narrow emission lines forming at the footpoints of accre- 
tion funnels (e.g. Ca II IRT, He I and Fe II emission lines) are 
mainly modulated by rotation. Broad emission lines (e.g., 
Balmer lines or the broad emission component of the IRT) 
are also modulated by rotation but apparently include a 
higher level of intrinsic variability than the narrow emis- 
sion lines. Spectropolarimetric observations spanning (and 
densely sampling) at least 2 complete rotation cycles are 
required to study in more details the relative strength of ro- 
tational modulation and intrinsic variability in most spectral 
lines of interest. 

Rotational modulation of photospheric and narrow 
emission lines can be mostly attributed to the presence 
of dark spots and hot accretion regions at the surface of 
BP Tau. We find that both types of features grossly overlap, 
concentrate near the pole and cover altogether a relative area 
of about 2%, similar to what was observed on V2129 Oph 
(D07). This agrees with previous independent observational 
and theoretical studies estimating that the accretion spots 
typically cover no more than a few percent of the stellar sur - 
face Ce.g.. lValenti fe Johns-Krullll2004 IJardine et aLll2006T l. 
It also agrees with the findings of lValenti et al.l (2003) that 
BP Tau (and other cTTSs) are faintest during times of max- 
imum emission fluxes and polarisation signatures. 

Veiling is also detected in the photospheric spectra, at 
an absolute level varying between and 60%, with relative 
changes of up to ±15% within a single run. We find that 
veiling only weakly correlates with the amount of emission 
and polarisation, e.g., in He I lines, and is subject to intrin- 
sic fluctuations much larger than those seen in line emission 
fluxes and po l arisat ion signatures. This is similar to what 
IValenti et all (2003) report, strengthening their conclusion 
that veiling is rather due to stochastic variations in accre- 
tion rates rather than to the magnetic geometry and the 
location of accretion spots. It also suggests that veiling (as 
a pure continuum excess) is likely not causing the observed 
Stokes / profile distortions in photospheric lines (as both 
would otherwise neatly correlate together), bringing further 
support to our model in which profile distortions of photo- 
spheric lines are mainly caused by dark surface spots. 

If our model is confirmed, it would indicate that that 
the heat produced in the accretion shock is not transferred 
to the photosphere efficiently enough to warm it up above 
the temperature of the non-accreting photosphere. 

6.2 Magnetic topology and origin of the field 

Stokes V signatures are detected on BP Tau, both from pho- 
tospheric and narrow emission lines; the very different Zee- 
man signatures that both sets of lines exhibit (featuring in 
particular opposite line-of-sight polarities) suggest that they 
trace magnetic fields from different complementary regions 
of the stellar surface, with emission lines tracing the fields 
of accreting regions and photospheric lines tracing the fields 
of non-accreting regions (as in V2129 Oph, see D07). With 
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Figure 15. Magnetospheric topology of BP Tau as derived from potential extrapolations of the Feb06 surface magnetic field distribution 
(top panel of Fig. 1140 . The magnetosphere is assumed to extend up to the inner disc radius, equal to 3.5 and 7.5 R+ in the left and 
right panels respectively. The complex magnetic topology close to the surface of the star is very obvious. In both cases, the star is shown 
at rotational phase 0.0. The colour patches at the surface of the star represent the radial component of the field (with red and blue 
corresponding to positive and negative polarities); open and closed field lines are shown in blue and white respectively. 
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Figure 16. Location of the footpoints of accreting field lines at the surface of the BP Tau (white circles), assuming that the magnetosphere 
(extrapolated from the Feb06 magnetic maps) extends to 7.5 R*, and that the accretion disc coincides with the stellar rotation equator 
(no magnetic warp, left panel) or with the stellar dipole magnetic equator (magnetic warp, right panel). The colours at the surface of 
the star depict the radial field component (with red and blue corresponding to positive and negative polarities). Note that phase runs 
backwards with longitude, i.e., decreases from 1 to while longitude increases from 0° to 360°. 



this model, we are able to recover the large-scale magnetic 
topology of BP Tau from both sets of Zeeman signatures, 
and at both observing epochs. The reconstructed magnetic 
field involves mainly a dominant dipole component of 1.2 kG 
and a strong octupolar term of 1.6 kG, both only slightly 
tilted with respect to the rotation axis. 

Zeeman signatures demonstrate that strong fields are 
present at the surface of BP Tau. We find that the aver- 
age magnetic filling factor, describing the relative amount 
of light from magnetic regions within each local surface 
area at the surface of BP Tau, is about 25%; in this con- 
text, the peak magnetic fluxes of 3 kG we reconstruct on 



BP Tau indicate that magnetic fields likely reach maximum 
field strengths of up to 12 kG, mostly in accreting regions. 
This result is supported by the strong longitudinal fields 
(in excess of 8 kG) traced with the He I 667.815 nm line, 
presumably forming mostly in the magnetic accreting re- 
gions (with little contribution from adjacent non-magnetic 
areas) and thus less prone to longitudinal dilution than other 
narrow emission lines. This is also in qualitative agreement 
with th e study of magnetic broad ening of infrared lines of 
BP Tau lJohns-Krull et al.l (|l999bl) . indicating that fields of 
up to 10 kG are likely to be present at the surface of the star. 
We note that the magnetic pressure within the accreting re- 
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gions vastly exceeds the gas pressure of the surrounding non- 
magnetic photosphere, making it difficult to understand why 
and how accreting spots are spatially confined (even more 
so than in the case of sunspots); one possible (though fairly 
speculative) reason could be that magnetic fields are also 
present in the quiet photospheric regions of BP Tau, but 
only in the form of small-scale tangled multipolar regions 
producing no detectable circular polarisation. 

Considering the whole surface of BP Tau, we find that 
the average magnetic flux over non-accreting regions is 
about 1.2 kG while that over the complete star is 1.4 kG; 
looking now at the visible hemisphere only (e.g., at phase 0.0 
in the Feb06 image) , we find that the average magnetic flux 
from non-accreting regions is again about 1.2 kG whereas 
that over the full visible hemisphere is 1.8 kG. The first 
of these 2 values is in good agreement with average mag- 
netic fl uxes derived from optical lines by I Johns-Krull et al.l 
(1999b). These authors however find that infrared lines in- 
dicate l arger average magn etic fluxes (of as much as 2.8 kG, 
though Ijonns-Krulll [20071 mentions 2.2 kG only), presum- 
ably because these lines also trace magnetic fields in cool 
highly magnetic spots whose brightness contrast with re- 
spect to the photosphere is much smaller at infrared than 
at optical wavelengths. Although we also find that strongest 
magnetic fields indeed concentrate in cool spots, the average 
brightness-unweighted magnetic flux we derive (1. 8 kG) is 
sma ller than that estim ated by Ijohns-Krull et all l|l999bl ) 
and I Johns- Krulll l|2007h . suggesting that we still likely miss 
magnetic flux from cool spots (e.g., from dark non-accreting 
small-scale magnetic regions). 

Comparing with V2129 Oph, we find that BP Tau hosts 
a 4 times stronger dipole field and a half as strong toroidal 
field. We speculate that this difference likely reflects the fact 
that BP Tau is still fully-convective, while V2129 Oph had 
recently started to build up a radiative core (D07) . Fully con- 
vective stars are indeed much more successful at triggering 
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in preparation) th an their more mas sive partly-convective 
counterparts (e.g.. iDonati et al"1l2003T l. If confirmed, this re- 
sult would suggest that the strong fields of low-mass cTTS 
are likely dynamo-generated rather than fossil fields; the un- 
derlying processes capable of producing such fields almost 
without the help of differential rotation a re not yet fully un- 
derstood from a theoretic al point of view l|Chabrier fc Kiikerl 
2006; Dobler et al.ll2006h . In particular, it is not clear how 
dynamo processes are able to produce magnetic fields that 
vastly exceed the thermal equipartition value, such as those 
of BP Tau and fully convective M dwarfs. 

The large-scale magnetic topology of BP Tau has ap- 
parently undergone no more than small changes between 
Feb06 and Dec06 - apart from a global phase shift of about 
0.25 rotation cycle that could be due to a slight error on 
the assumed rotation period. It suggests that the lifetime 
of BP Tau's magnetic topology is comparable to the times- 
pan between our two sets of observations, i.e., much longer 
than those of partly-convective active stars (whose field gen- 
erally chan ges beyond recogni tion in no more than a few 
weeks, e.g.. |Ponati et al.ll2003f ). This result is in agreement 
with recent claims that magnetic fields of fully-convective 
main-sequence dwarfs are stable on timescales of at least 



1 yr (|Morin et alj|2007l . Morin et al., 2008 in preparation). 
If confirmed, it would imply that the nominal rotation pe- 
riod of BP Tau (7.6 d) is slightly underestimated; a better 
phase match between both images (separated by 298 d or 39 
rotation cycles) is obtained when assuming that the rotation 
period is 7.65 d. Note that this new esti mate is very close 
to the value of 7.64 d initially derived bv lVrba etai] (|l986l ) 
from multicolour photometric data. 

6.3 Disc-star magnetic coupling 

Several theoretical papers (e.g. IKonigll Il99ll; IShu et al.l 
1 19941 ; ICameron fc Campbelllll993l ; lLong et al.ll2005l ) studied 
how the stellar magnetic field interacts with the surrounding 
accretion disc and disrupts its vertical structure close to the 
star. They showed further that the balance between accre- 
tion torques and angular momentum losses causes the rota- 
tion of the star to evolve towards an equilibrium in which 
the disc disruption radius lies close to ta and just inside 
the co-rotation radius re- They proposed that this coupling 
causes cTTSs to slow down to the Keplerian orbital period 
at a radius about 10-50% larger than r\ (i.e., r\ — 0.8 rc), 
explaining why cTTSs are on average rotating more slowly 
than their disc-less equivalents. This scenario is refered to 
as 'disc-locking' in the literature. 

For BP Tau, we obtain that the 1.2 kG large-scale dipole 
field we observe yields an equilibrium radius rc ranging be- 
tween 6.3 and 8.9 R* (depending on whether magnetic dif- 
fusivity is due to buoyancy or turbulent diffu sion) according 
to the model of lCameron fc Campbelll (1 19931 ) . Simi lar results 
are obtained from the model of Long et al.l (|2005l ), yielding 
rc = 8.5 R t and ta — 6 ii* when applied to BP Tau. This is 
in agreement with both our own independent estimates of rc 
(about 7.5 R *) and t a (lar ger than 4 see Sec. [5}. Using 
the model of lK5nigl (199l|) however, we find an equilibrium 
radius of 2.7 in strong contradiction with both our es- 
timates of rc and ta- As for V2129 Oph, we find that the 
disc-locking scenario is compatible with observations when 
the magnetic coupling between the st ar and its accretion disc 
is des crib ed by models like those of ICameron fc Campbelll 
i| 19931 ) or lLong et all l|2005h . 

Our results suggest that winds on BP Tau are obvi- 
ously not strong enough to blow open field l ines larger than 
3 i?* (e.g.. lSafierlll998l ; lMatt fc Pudritzll2004l ) and to prevent 
disc magnetic disruption to take place and magnetospheric 
accretion to occur as far as ta- If it were the case, accre- 
tion spots would distribute differently over the surface of 
the star, with a significant fraction of mass being accreted 
towards the equator, in strong contradiction with our ob- 
servations. Our observations also suggest that the accretion 
disc is magnetically warped at the vicinity of ta\ a similar 
conclusion was reached in the c ase of the cTTS AA Tau on 
completely different arguments (|Bouvier et al.|[2007h . 



7 CONCLUSION 

In this paper, we report the detection of spectropolarimet- 
ric Zeeman signatures on the cTTS BP Tau using mostly 
ESPaDOnS /CFHT (with help from NARVAL/TBL). Cir- 
cular polarisation signatures in photospheric lines and in 
narrow emission lines tracing magnetospheric accretion are 
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monitored throughout most of the rotation cycle of BP Tau 
at two different epochs (Feb06 and Dec06). We find that 
polarised and unpolarised spectral proxies tracing the pho- 
tosphere and the footpoints of accretion funnels show tem- 
poral variations that are mostly attributable to rotational 
modulation. 

From our spectropolarimetric data sets of photospheric 
and narrow Ca n emission lines simultaneously, we recon- 
struct, using tomographic imaging, the large-scale magnetic 
topology and the location of accretion spots at the surface of 
BP Tau at both epochs. We find that the magnetic topology 
of BP Tau involves a dominant (1.2 kG) dipole but also a 
strong (1.6 kG) octupole, both slightly (though differently) 
tilted with respect to the rotation axis. In particular, the 
strong dipole component makes BP Tau fairly different from 
V2129 Oph, another (more massive) cTTS on which a sim- 
ilar study was carried out (D07). Accretion spots coincide 
with the two main high-latitude octupole poles and overlap 
with dark photospheric spots; they each cover about 2% of 
the stellar surface. 

Despite clear variability between both epochs, the large- 
scale magnetic topologies we reconstruct are nevertheless 
grossly similar, suggesting an overall lifetime longer than 
6 months. The strong mostly-poloidal, nearly-axisymmetric 
field of BP Tau (and in particular its long-lived dipole com- 
ponent) is very re miniscent of magnet i c topologies of fully- 
convective dwarfs llDonati et a l. 2006a; Hal linan et al.ll2006l ; 
iMorin et ai1l2007l ; lHallinan et al.ll2007l . Morin et al., 2008 in 
preparation). It suggests that the strong large-scale poloidal 
fields hosted by fully-convective cTTSs such as BP Tau 
(but absent in more massive non fully-convective cTTSs like 
V2129 Oph) are likely not fossil remants, but rather result 
from vigorous dynamo action operating within the bulk of 
their convective zones. 

Preliminary modelling suggests that the magnetosphere 
of BP Tau must extend to distances of at least 4 7?* to pro- 
duce accretion spots at a latitude roughly matching those 
we observe. At the very least, it demonstrates that mag- 
netic field lines from the protostar are not blown open close 
to the surface by a stellar wind, but are apparently capa- 
ble of coupling to the accretion disc beyond 3 ii*. Our esti- 
mates of the magnetospheric and coro tation radii for BP Tau 
match the theoretical predict ions of lCameron fe Campbell! 
1 19931 ) and lLong et al.l l)2005r ). thus strengthening the idea 
that star/disc magnetic coupling may indeed be responsi- 
ble for the slow rotation of fully-convective cTTSs such as 
BP Tau. 
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